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Abstract 

The purpose of this paper is to analyze how the disorder affects the dynamics 
of critical fluctuations for two different types of interacting particle system: 
the Curie-Weiss and Kuramoto model. The models under consideration are 
a collection of spins and rotators respectively. They both are subject to a 
mean field interaction and embedded in a site-dependent, i.i.d. random envi- 
ronment. As the number of particles goes to infinity their limiting dynamics 
become deterministic and exhibit phase transition. The main result concern 
the fluctuations around this deterministic limit at the critical point in the 
thermodynamic limit. From a qualitative point of view, it indicates that when 
disorder is added spin and rotator systems belong to two different classes of 
universality, which is not the case for the homogeneous models (i.e., without 
disorder). 
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1 Introduction 

Interacting particle systems with mean field interaction are characterized by the 
complete absence of geometry in the space of configurations, in the sense that the 
strength of the interaction between particles is independent of their mutual position. 
The advantage of dealing with this kind of models is that they usually are analyt- 
ically tractable and it is rather simple derive their macroscopic equations. Even 
if the mean field hypothesis may seem too simplistic to describe physical systems, 
where geometry and short-range interactions are involved, mean field models have 
been recently applied to social sciences and finance, as in [3 [101 O [15] . 
We briefly introduce the general framework and some of its peculiar features. By 
mean field stochastic process we mean a family x^ N ' — (x^ N ' (t))t>o with the fol- 
lowing characteristics: 

• x( N \t) = (x[ N \t), X2*\t), . . . , x^\t)j is a Markov process with N compo- 
nents, taking values on a given measurable space (E,£); 

• Consider the empirical measure 

1 N 

:= ^E^< N) (*)' 
fc=i 

which is a random probability on (E, £). Then (pN(t))t>o is a measure- valued 
Markov process. 

Although this is by no means a standard definition of mean field model, it captures 
the basic features of the specific models we will consider. 

Let (F, F) be a topological vector space, and h : E — > F be a measurable function. 
Objects of the form 

r 1 N 

hd PN (t) = -j2h(4 N) (t)) 

J fc=l 

are called empirical averages. In the case the flow (J hdptq{t)) t >Q is a Markov 
process, we say J hdp^{t) is an order parameter. Note that the empirical measure 
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itself is an order parameter (taking F = set of signed measures on E, and h(x) = 5 X ). 
Whenever possible, it is interesting to find finite dimensional order parameters, i.e. 
order parameters for which F is finite dimensional. 

One of the nice aspects of mean field models is that, in many interesting cases, one 
can prove a Law of Large Numbers (as N — > +00) for the order parameters, and 
characterize the deterministic limit as a solution of an ordinary differential equation. 
This limit is often called the McKean- Vlasov limit. In particular, the differential 
equation describing the limit evolution of the empirical measure, will be referred to 
as the McKean- Vlasov equation. This equation has the form 

d r 

di q = ^ 

where £ is a nonlinear operator acting on signed measures on E (even though other 
spaces may be more convenient for the analysis of C). 

Our main interest is the study of the fluctuations of the order parameter around 
its limiting dynamics. We can capture different features of these fluctuations de- 
pending on whether or not the time is rescaled with N. If time is not rescaled 
and we consider the evolution in a time interval [0,T], with T fixed, a Central 
Limit Theorem holds for the order parameter for all regimes; in other words, the 
fluctuations of the order parameter converge to a Gaussian process, which is the 
unique solution of a linear diffusion equation. Whenever time is rescaled in such 
a way T goes to infinity as N does, we may observe different behaviors. To avoid 
further complications, we assume the Markov process x^ N \t) has a "nice" chaotic 
initial condition: x^\o) , x^\o) , . . . ,x^\o) are i.i.d. with common law qo(dx), 
where qo is a stationary, locally stable solution of the McKean- Vlasov equation (the 
system is in local equilibrium) . 

• Subcritical regime. Suppose qo is the unique stationary solution of the McKean- 
Vlasov equation, and it is linearly stable (i.e. stable for the linearized equa- 
tion). Then we expect the Central Limit Theorem holds uniformly in time; 
in particular, this provides a Central Limit Theorem for the stationary distri- 
bution of Some results in this direction are shown in [T2] . 

• Supercritical regime. Suppose the set of stationary, linearly stable solutions 
of the McKean- Vlasov equation has cardinality greater than 1. In this case 
metastability phenomena occur at a time scale exponentially growing in N. 

• Critical regime. This is the case in the boundary of the subcritical regime: 
denoting by £ the linearization of C around goj the spectrum Spec(£) of £ is 
contained in {z E C : Re{z) < 0}, but there are elements on Spec(£) with zero 
real part. Under a suitable time speed-up, the elements of the corresponding 
eigenspaces may exhibit large and, possibly, non-normal fluctuations (see [S] 

m- 

Of course the three regimes described above do not cover in general all possibilities, 
since stable periodic orbits or even stranger attractors may arise. Moreover, the 
same model could be in different regimes depending on the values of some parame- 
ters (phase transition). 

The main subject of this paper is the analysis of the dynamics of the critical fluc- 
tuations in disordered mean field models. 

We consider a mean field model and we add a site-dependent, i.i.d. random environ- 
ment, acting as an inhomogeneity in the structure of the system; we aim at analyzing 
the effect of the disorder in the dynamics of critical fluctuations, as compared with 
the homogeneous case. We deal with the Curie- Weiss and the Kuramoto models. 
We are not aware of similar results concerning non-equilibrium critical fluctuations 
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in presence of disorder. Static fluctuations for the random Curie- Weiss model have 
been studied in pQ. 

We now give the basic ideas of how the dynamics of critical fluctuations are deter- 
mined. As we mentioned above, the deterministic limiting dynamics of the order 
parameter is described by a nonlinear evolution operator C. The linearization of 
this equation around a stationary solution gives rise to the so called linearized op- 
erator £. This operator is also related to the normal fluctuation of the process. At 
the critical point this operator has an eigenvalue with zero real part, while all other 
elements of the spectrum have negative real part. The eigenspace of the eigenvalue 
with zero real part will be called critical direction, and usually happens to have low 
dimension: critical phenomena involve the empirical averages corresponding to this 
subspacc. Thus, our analysis follow the following points. 

• Locating the critical direction. 

• Determining the correct space-time scaling for the critical fluctuations. This 
requires an approximation of the time evolution of the order parameter that 
goes beyond the normal approximation. 

• Proving that the rescaled fluctuations vanish along non-critical directions. 
This will be done using the method of "collapsing processes" : it was developed 
by Comets and Eisele in [5] for a geometric long-range interacting spin system 
and was previously applied to a homogeneous mean field spin-flip system in 

HE- 

• Determining the limiting dynamics in the critical direction. It will be done 
using arguments of perturbation theory for Markov processes, which has been 
treated in [17] . and of tightness, applied to a suitable martingale problem. 

From a qualitative point of view, our results indicate that when disorder is added, 
spin systems and rotators belong to two different classes of universality, which is 
not the case for homogeneous systems. Roughly speaking, in spin systems the 
fluctuations produced by the disorder always prevail in the critical regime: these 
fluctuations evolve in a time scale of order N~i, while the critical slowing down for 
homogeneous systems is N?. For rotators, the disorder does not modify the 
slowing down. However, as the "strength" of the disorder increases, the Kuramoto 
model undergoes a further phase transition: for sufficiently small disorder, the 
dynamics of critical fluctuations converge to a nonlinear, ergodic diffusion, as in the 
homogeneous case; for larger disorder, the limiting diffusion loses ergodicity, and 
actually explodes in finite time. 

We finally remark that in [?] we have analyzed the critical fluctuations for a spin 
system close in spirit to the Curie- Weiss model, although with a less general disorder 
distribution. 

2 The Random Curie- Weiss Model 
2.1 Description of the Model 

Let 5? = { — 1,+f} be the spin space, and \i be an even probability on HL Let 
also r\ = & R N be a sequence of independent, identically distributed random 

variables, denned on some probability space (fi, J 7 , P), and distributed according to 
/i. They represent a random, inhomogeneous magnetic field. 

Given a configuration a = {(Jj)f =1 € S? N and a realization of the magnetic field 77, 
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we define the Hamiltonian H N (a,rj) : y N x R N -> R as 

#w (2, 2) = - ^ -PzLvjVj, (!) 
j,fe=i j=i 

where <7j is the spin value at site j, and 77^- is the local magnetic field associated with 
the same site. Let j3 > be the inverse temperature. For given r\, a(t) — ((Jj(t))^ =1 , 
with t > 0, is a TV-spin system evolving as a continuous time Markov chain on 5? , 
with infinitesimal generator Lm acting on functions / : — > R as follows: 

N 

LnKct) = XI e-^'^+^VJ/Qr), (2) 

where VJ/(ct) = /(aj ? ) — / (a;) and the fc-th component of gj , which is the spin flip 
at the site j, is 

a i _ f o-fe for fc ^ j 
k 1 — Ofc for k = j 

The quantity e _/3o ' J ( m « + ' ?J * 1 represents the jump rate of the spins, i.e. the rate at 
which the transition Oj — > —<jj occurs for some j. The expressions ([T]) and 
describe a system of mean field ferromagnetically coupled spins, each with its own 
random magnetic field and subject to Glauber dynamics. The two terms in the 
Hamiltonian have different effects: the first one tends to align the spins, while the 
second one tends to point each of them in the direction of its local field. 

Remark 2.1. For every value of 77, @ has a reversible stationary distribution pro- 
portional to exp[— Hn(o_, 77)]. 

For simplicity, the initial condition <r(0) is such that (<Xj(0), r]j)^ =l are independent 
and identically distributed with law A. Note that, since the marginal law of the rjj's 
is /i, A must be of the form 

A(cr, drj) = q (cr, T])n(dri) (3) 

with go(l, v) + <7o(~ 1, v) = 1> /J-almost surely. The quantity (uj (t)) tg [ T j represents 
the time evolution on [0, T] of j-th spin value; it is the trajectory of the single ji-th 
spin in time. The space of all these paths is D[0,T], which is the space of the 
right-continuous, piecewise-constant functions from [0,T] to 5? . We endow 2?[0,T] 
with the Skorohod topology, which provides a metric and a Borel c-field (see [TT] 
for details). 



2.2 Limiting Dynamics 

We now describe the dynamics of the process (J2J), in the limit as N — > +00, in a 
fixed time interval [0,T]. Later, the equilibrium of the limiting dynamics will be 
studied. These results are special cases of what shown in [6], so proofs are omitted. 
More details can also be found in (3J. 

Let (aj[0,T])f =1 G (D[0,T]) N denote a path of the system in the time interval [0,T], 
with T positive and fixed. If / : 5? x R — > M, we are interested in the asymptotic 
(as N — > +00) behavior of empirical averages of the form 

N 
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where (/9/v(^))te[o,Tl is the flow of empirical measures 

N 

X 



PN{t) ■= 



We may think of pn '■= (Aiv(i))te[o,T] as a cadlag function taking values in M.\{Sf x 
R) , the space of probability measures on 5? x R endowed with the weak convergence 
topology, and the related Prokhorov metric, that we denote by dp{ ■ , • ). 
The first result we state concerns the dynamics of the flow of empirical measures. 
We need some more notations. For a given q : x M — > R, we introduce the linear 
operator C q , acting on / : J? x K — > R as follows: 

£,/(<r,»7) := V ff [e-^( ra «^/((r,»7) 
where 

m 9 := J [q(l, rj) - q(-l, r?)] /i(c^). 

Given 77 e R w , we denote by 7^ the distribution on (PfOjT])^ of the Markov 
process with generator @ and initial distribution A. We also denote by 

V N {da[0, T],drj) := V% (da[0, T]) fi® N (*/) 

the joint law of the process and the field. 

Theorem 2.1. 27ie nonlinear McKean-Vlasov equation 

qo(<r,r)) given in (J3J) 

admits a unique solution in C 1 [0, T], , arte? qt(',v) * s probability on 5? , 

for p-almost every r\ and every t > 0. Moreover, for every e > £/iere exists 
C(e) > smc/i £/ia£ 

^ I sup d P ( PN {t),q t ) >e ) <e~ c ^ N 
\te[o,T] J 

for N sufficiently large, where, by abuse of notations, we identify q t with the prob- 
ability <7t(cr, rj)fj,(dn) on x R. 

Thus, equation Q describes the infinite-volume dynamics of the system. Next 
result gives a characterization of stationary solutions of (|4]). 

Lemma 2.1. Let q* : y x I -> 1, such that q*(a, •) is measurable and q*{-,i]) is 
a probability on 5f . Then <?* is a stationary solution of i.e. C q ,q* = 0, if and 
only if it is of the form 

g * (a ^ 2cosh(/3( TO , +?? )r (5) 
where satisfies the self-consistency relation 

[g*(l,?7)-g*(-l,»7)]/x(d»7). (6) 



Moreover, = is always a solution of ([6]) and it is linearly (resp. neutrally) 
stable if and only if 

P [ -J*™ < (resp. =) f. (7) 
J cosh [f3r]) 
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Remark 2.2. The transition between uniqueness and non- uniqueness of the solution 
of ([6]) in general is not related to the change of stability for to* = 0. If the distri- 
bution p, is unimodal on R, the two thresholds coincide: the paramagnetic solution 
is linearly stable when it is unique and unstable when it is not. In case we choose 
M = i (<^7 + ^— w hh V > 0j the phase diagram is more complex: when ([7J fails, 
the paramagnetic solution of (jSJ) is either unstable, and it coexists with a pair of op- 
posite stable ferromagnetic solutions, or may recover linear stability, coexisting with 
a pair of unstable ferromagnetic solutions and a pair of stable ferromagnetic ones 
(see [5] for details). A more general p may give rise to arbitrarily many solutions 
of ©. 



2.3 Dynamics of Critical Fluctuations I j3 



cosh 2 (firf) 



The results of this section are concerned with the fluctuation flow 

p N (t) :=y/N[p N (t)-q t }, 



(8) 



that takes values on the space of signed measures on Yxl. It is very convenient 
to assume that the process starts in local equilibrium, i.e. qo(o~, 77) = q*{o~, 77), where 
(a, 77) is a stationary solution of Q ; it should be not hard to extend all next 
results to a general initial condition. The proofs of all results stated here will be 
given in Section [5J We first state results valid for all temperatures; later, Lemma 
2.31 Proposition 12.31 Theorems 12.21 and 12.31 are restricted to the critical case. 
Functions from / xl are all of the form F(a, rj) = 7(77) + o<j>{rj). However 

N 



j(rj)dp N (t) = VN 



3=1 



does not change in time, and has a Gaussian limit for every 7 £ L 2 (/i). Thus, we 
are only interested in the evolution of integrals of the type 



a(f>(r))dp N (t). 



It is therefore natural to control the action of the generator Ln on functions of a 
and rj of the form ip (f<T(f>(jf)dptf), with 



Pn ■= 

Proposition 2.1. Let ip : K 

the measure on K defined by 



1 N 



Vj) 



3=1 



v(dri) 



be of class C 1 , and cf> € (L 2 (i/))™ ; where v is 
H(dq) 



cosh(/3(TO* + 77)) 



(9) 



Ther. 



L N ^ I J <J(j)(r))dp N 
= 2j2dii>([ o-<f>(r,)dp 

+ 2t 8 J«(/*)*)/ c ^ K+ii)) 

* 1 .7 — 1 



sinh(/3(TO* + rf))4>i{ri)dp N - J aZ^r^dpN 

„(*,)+ o(l), (10) 
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where 

= cosh(/3(m, + V )W,V) ~P I ,,t^ , ■ (11) 



cosh(/3(m» + ??)) 
Moreover the remainder o(l) m 175)) is 0/ f/ie /orm 



R N U H(a,r))dp N ) (12) 
where H(a,rf) is the vector-valued function 

H{a,rj) = {a<j}{rf),(T, [cosh(/3(m» +17)) - a-sinh(/J(m» +?7))]^(r?), 

[<7COsh(/3(m» + 77)) - smh(/3(m* + v))]^)) > 



lim sup Rn(x,u,z,w) = (13) 



/or every M > 0. 



Proposition 12. H whose proof consists of a rather standard computation that will be 
sketched in Section^ is the essential ingredient for proving a Central Limit Theorem 
for the empirical flow, i.e. to show that the fluctuation flow converges in law to 
a Gaussian process. The proof of this result requires to identify an appropriate 
Hilbert space for the fluctuations pjy (see e.g. [5] for related results). Our main aim 
is, however, to describe large-time fluctuations at the critical points; the additional 
technical difficulties arising, have not allowed us to obtained the desired results 
under the present assumptions, in particular with no requirements on the field 
distribution p (except evenness). Thus we find it preferable to make the following 
assumption at this point. 

(F) fi has finite support @. 

Under assumption (F), the space L 2 {v) is finite-dimensional. Together with the 
following simple result, this greatly simplifies the analysis of fluctuations. 

Lemma 2.2. The operator £ defined in (|11[) is self-adjoint in L 2 (v). 

Now, for m := |supp(/i)|, let ipo, ipi, . . . , ip m — \ be a complete set of eigenvectors 
for £, with eigenvalues Ao < Ai < . . . < A m _i. Proposition 12.11 together with 
the classical Corollary 8.7, in Chapter 4 of [TT], yields the following Central Limit 
Theorem, whose standard proof is omitted. 

Proposition 2.2. Setxf\t):= J aip l (rj)dp N (t) . Then, underV N , (x^)™ 

V / i— 

converges in law to the Gaussian process (Xi)™^ 1 solving the following linear stochas- 
tic differential equations 

dX t {t) = [Mi - XiXi(t)] dt + h dWi(t) 

where 

• (Xq(0), Xi(0), . . . ,X m _i(0), J%o, J4f[, . . . , J%n-i) is a centered Gaussian vector 
with 



Cov 



(x l (o),x J (o)) = J ipMwWvW) 

-J tanh( / 9(m* + r]))fj,(dn) J fi(v) tanh(/3(m* + ri))n(dri) 
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Cov 



Cov 



= f<Pi(v)Vj(v) smh 2 (/3(m» + rf))fi(dri) 

-J <fii(v) sinh03(m* + n))fi(dn) J <Pj(ri) smh(/3(m, + n))fl(dn) 

(J%,Xj(0)) = J <Pi{ri)<pj{ri) sinh(/3(m, + r,)) tanh(/3(m» + rj))n(dr)) 

-J Viiv) sinh(/3(m» + r}))fi(dr}) J <fij(ri) tanh(/9(m* + r)))fl(dn) 

• (Wi)™^ 1 are independent standard Brownian motions, that are independent of the 
vector pf (0),Xi(0), . . . , X m _i(0), J#o\ • • • 

Note that the randomness of the field persists in the limiting dynamics of fluctu- 
ations, due to the correlated, constant random drifts M'i. Observe that M[ = if 
/.t = (So, i-e. when the random field is absent. 

We now look more closely at fluctuations around the paramagnetic solution m* = 
at the critical regime, i.e. for those values of /3 for which (3 CO £$'(j3ri) = 1. 

Lemma 2.3. Assume f3 = 1 and, to* = 0. Then £ is nonnegative, and 

its kernel is spanned by the function g^ggrg^ ■ 

In the critical regime (3 ^^npfi = 1, we have Ao = 0, and \i > for i > (it 
is actually easily shown that A; > 1 for i > 0). It follows that the process Xo(t) in 
Proposition 12.21 has a variance that diverges as t — > +oo. A sharper description of 
the large time fluctuations is obtained by considering more "moderate" fluctuations: 

p N : = N~^p N . 

The following result improves the expansion given in Proposition 12. II 
Proposition 2.3. Under the same assumptions of Provosition lKTi and the further 
conditions (3 C mifi(j)rj) = 1 cmd to* = 0, we have 

L N i> U o-<j>(ri)dp N \ = L (1 V 

+ 2JV-*5^0i^ (J a^dpN^j J smhd3r})<t)i{ri)dp 



N 



+ N~^L^ip + N~^L^ip + o (n~*) , (14) 



where 



L^ip := ~2j2d.^ o-<j>{rj)dp N ^j J a£Mv)df>N 



-2/3^ diip (^J acf)(ri)dpNj J adp N J a smh(Prj)4>i{ri)dp 
£ (3) V> := y^dji/j {^J o-^rftdpN^j 2/3 J cosh{l3n)(j) t (n)dp N J adp N 



d2 



2 f ft 3 f 4>i(v) ( f x 3 

adp N ) J a cosh(l3ri)(t>i{rf)dp N + — J — -y /^(<iT7) I J adp N 



2 E (/ ^fow*) / 



cosh(/3?7) 
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Moreover the remainder o a j in Jj^P * s of the form N 2 wi£/i iijy satisfying 



Note that in Proposition l2.31 functions depending only on 77 are still integrated with 
respect to p, rather than p; indeed, by the standard Central Limit Theorem, those 
integrals with respect to p have a Gaussian limit under Vn- 

Proposition 12.31 allows to deal easily with the homogeneous case p = Sq. Using the 
notations of Proposition 12.21 we have m = 1, ipo = 1. Thus, using Proposition 
with n = 1, <f> = 1 and (3 = f3 c = 1, we easily observe that L^ijj = L^tp = 0, and 



i (3 V = -\ (J odpN^J V (J o-dp N \ + 2^ (J adp N 



Using convergence of generators as in Proposition ^ . 2l we readily obtain the dynamics 
of large-time critical fluctuations for the homogeneous model. This result is a simple 
special case of what obtained in [5]. 

Theorem 2.2. Assume p = Sq, and f3 = 1. The stochastic process 

Y N (t) := fadp N (VNt) 



converges weakly, under Vn, to the unique solution of the stochastic differential 
equation 

dY(t) = -lY 3 (t)dt + 2dW(t) 
Y(0) = 

where W is a standard Brownian motion. 

As we will see (proofs are in Section [5]), the inhomogeneous case requires more 
sophisticated arguments. 

Definition 2.1. We say that a sequence of stochastic processes (£n(t)) n j for t G 
[0, T], collapses to zero if for every e > 0, 

lim P\ sup |f n (t)| >e) =0 

y te [o,r] y 

Theorem 2.3. Assume — 0, (3 ^ffijs^ — U an d, for i — 0, 1, . . . , m — 1, let 

Y^ N \t) := f a 9i {n)d~PN{Nhl (15) 



where ipo, . . . , <y9 m _i is the basis introduced in Proposition^^ Under Vn the pro- 



((N) \ rn 1 ( N) 

Yj (t)J collapse to zero, while Y" (tj converges in law to the process 

Y (t) := 2,m, 

where is a Gaussian random variable, with zero mean and variance tanh 2 (f3r]) p{drj) . 

Thus, the disorder has a dramatic impact on fluctuations at the critical points: 
fluctuations arise at a much shorter time scale (N* rather that N?), and have the 
simple form of a linear function with random slope. 
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3 The Random Kuramoto Model 



3.1 Description of the Model 

Let / = [0, 27r) be the one dimensional torus, and /i be an even probability on R. Let 
also T] = (r]j)jLi £ K N be a sequence of independent, identically distributed random 
variables, denned on some probability space (fi,.F, P), and distributed according 
to /i. Given a configuration x — (xj)^ =l £ I and a realization of the random 
environment ?y, we can define the Hamiltonian H^(x,rj) : I x R N -> K as 

Q N N 
j,k=l J=l 

where Xj is the position of the rotator at site j and u>r)j, with uj > 0, can be 
interpreted as its own frequency. Let 9, positive parameter, be the coupling strength. 
For given r), the stochastic process x(t) — {xj(t))^ =ll with t > 0, is a iV-rotator 
system evolving as a Markov diffusion process on I N , with infinitesimal generator 
Ln acting on C 2 functions / : I N — > R as follows: 

N o2 - ^ 




(17) 



1 N 

^ e,; *" = ^E e " J ' ( 18 ) 

where < < 1 measures the phase coherence of the rotators and \&jv measures 
the average phase. We can reformulate the expression of the infinitesimal generator 
1(17)1 in terms of (ITS)) : 

1 N d 2 f N df 

L N f(x) = - ^fe) + E ^ + 9rN s[n ^ N - X M ^fe) • ( 19 ) 

j=l 3 j=l ■? 



The expressions (1161) and ((19)) describe a system of mean field coupled rotators, 
each with its own frequency and subject to diffusive dynamics. The two terms 
in the Hamiltonian have different effects: the first one tends to synchronize the ro- 
tators, while the second one tends to make each of them rotate at its own frequency. 

For simplicity, the initial condition x(0) is such that (xj(0), f?j)jLi are independent 
and identically distributed with law A. We assume A is of the form 

\{dx, drj) = qo(x,7j)fi(dr])dx (20) 

with J j qo(x,r]) dx — 1, /i- almost surely. The quantity Xj(t) represents the time 
evolution on [0, T] of j-th rotator; it is the trajectory of the single j-th rotator in 
time. The space of all these paths is C[0,T], which is the space of the continuous 
function from [0, T] to /, endowed with the uniform topology. 
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3.2 Limiting Dynamics 

We now describe the dynamics of the process (fTTl) . in the limit as N — > +00, in a 
fixed time interval [0,T]. Later, the equilibrium of the limiting dynamics will be 
studied. These results are special cases of what shown in [5], so proofs are omitted. 
Let (xj [0, T])f =1 e (C[0, T]) N denote a path of the system in the time interval [0, T], 
with T positive and fixed. If / : I x K — > R, we are interested in the asymptotic (as 
N — » +00) behavior of empirical averages of the form 



N 

jyX;/fo(t),»fe)=: J fdPN® 



3=1 

where (pjv(^))te[o,Tl is the flow of empirical measures 

1 N 

PN(t) := J-fJ2 S {x 3 (t), Vj ) ■ 
3=1 

We may think of pjy := (pat (*))te[o,ri as a continuous function taking values in 
x K), the space of probability measures on / x M endowed with the weak 
convergence topology, and the related Prokhorov metric, that we denote by dp( ■ , • ). 
The first result we state concerns the dynamics of the flow of empirical measures. 
We need some more notations. For a given q : I x R — >• R, we introduce the linear 
operator C q , acting on / : / x R — > R as follows: 

Ar/M = ^(^^-^{H + ^sin^-x)]/^,^)}, (21) 

where 

rge 4 * 9 := J J e %x q(x , i]) fi(dn) dx . 

Given 77 e R^, we denote by T 3 ^ the distribution on (C[0,T]) N of the Markov 
process with generator (ITTl) and initial distribution A. We also denote by 

V N {dx[0, T},drj) := V% (dx[0, T]) fi® N (dry) 

the joint law of the process and the environment. 

Theorem 3.1. The nonlinear McKean-Vlasov equation 

dt ~ '-qtltyx, l\) 
qo(x,rj) given in (l20l) 



admits a unique solution in C 1 [[0, T], ^(dx (g> fj,)] , and qt{-,rj) is probability on I, 
for fi-almost every r\ and every t > 0. Moreover, for every e > there exists 
C{e) > such that 

V N ( sup d P ( PN (t),q t )>e) <e- c ^ N 
\te[o,T] J 

for N sufficiently large, where, by abuse of notations, we identify q t with the prob- 
ability qt(x, rj)ii(drj)dx on I x R. 

Thus, equation (|22[) describes the infinite-volume dynamics of the system. Since fj, 
is symmetric and the operator £ preserves evenness, we can suppose the average 
phase 1 iq t = 0, without loss of generality. Next result gives a characterization of 
stationary solutions of l[2"2"j). 
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Lemma 3.1. Let q* : I x R — > R 7 such that q*(x, •) is measurable and q*(-,rj) is a 
probability on I. Then q* is a stationary solution of (|22p . i.e. £ q „q* = 0, if and 
only if it is of the form 



+(i 



-2(ojT}y-\-0r* cosy 



dy 



, (23) 



where Z* is a normalizing factor and satisfies the self- consistency relation 



r* 



e M q*(x, rj) p(dn) dx . 



Moreover, r* = is always a solution of (|24l) and, letting 



1 + 4(w?7) 5 



(24) 



(25) 



we /iawe that 



1. if fi is unimodal on R, t/ien f/ie solution = is linearly (resp. neutrally) 
stable if and only if 9 < (resp. =) 9 C ; 

2. if fx — \(o~i + then the solution = is linearly (resp. neutrally) stable 
if and only if 9 < (resp. =) 9 C A 2 . 

Remark 3.1. The transitions uniqueness/non- uniqueness of the solution of (1241) and 
stability/instability of =0 in general do not occur at the same threshold. It does, 
however, in the case [T] of the previous Lemma. The phase diagram related to the 
case [2] is more complicated. We refer to [6] for further details. 

Remark 3.2. If r* = the stationary solution (|23|) reduces to q*(x,rj) := 



3.3 Dynamics of Critical Fluctuations 9 C 



li(dr]) 



1 + A(ur]) 2 _ 



-r 



The results of this section are concerned with the fluctuation flow 

p N (t) :=VN\p N (t)-qt], 



(26) 



that takes values on the space of signed measures on / x 1. It is very conve- 
nient to assume that the process starts in the particular local equilibrium qo(x, rf) — 
q*(x, if) = which is the stationary solution of (1221) corresponding to = 0. The 
proof of the Central Limit Theorem fProposition l3.2l) should be not hard also when 
q*(x,ri) is a sincronous stationary solution of (|22p . i.e. with =/= 0. The proofs of 
all results stated here will be given in Section [5] 

If is a function from / x R, we are interested in the evolution of integrals of the 
type 

4>(x,r])dp N (t) 



It is therefore natural to control the action of the generator Ljy on functions of x 
and 77 of the form ip (J <p(x, rj)dpiy) , with 



p N := ViV 



1 N 



= Vn 



1 N 1 
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Proposition 3.1. Let ip : R n ->• M &e 0/ class C 2 , and e (C 2 ([0, 2tt) x {-1, l})) r 
be 2ti -periodic in the first argument. Then 



L N 4> (^J (fiix^rfdpN^J = ^hdjip 



<j>(x,rj)dp N 



N% J dx 



£4>i(x,<n)dp N 

(x, n) sm(y - x)dp N dp N 



i,k=l 



where the operator 



d<t>i, ,d(j) k 



1 f dfa d(j) k 



(27) 



is the linearization of C, given by (|21l) . around the equilibrium distribution 



cos a; / cosy <f>(y, r/) dq* 

+ s'mx / sin y 4>(y,7]) dq* 



(28) 



Unlike the proof of Proposition 12.11 which requires an expansion of the generator, 
Proposition 13. II follows by the direct application of the generator; its proof is omit- 
ted. It provides the key computation for the proof of the Central Limit Theorem 
(Proposition 13.21 below). In order to simplify the analysis, we make the following 
assumption on the distribution of the random environment. 

(HI) n=±(6i + 6-i) 

Because of the structure of the system, it is reasonable to focus on functions from 
I x R of the forms <j)(x,rj) — cos(hx), sm(hx), r]cos(hx) or r/sm(hx), for ft, > 1 
integer, and thus on the behavior of 



X { ^' N) {t) := / cos(hx)dp N (t), X^' N) (t) := / sm(hx)dp N (t), 



xj^' N '(t) := / ncos(hx)dpN(t) and Xj^' N \t) := / r/sin(hx)dpN(t) 



Proposition ^. 1[ together with the classical Corollary 8.7, in Chapter 4 of [IT], yields 
the following Central Limit Theorem. 

Proposition 3.2. Assume (HI) holds. For r > 1, consider the following space of 
sequences 



H-r — ^ X — [ x h\ x h\ x h\ x \?^j ^ '■ ll x l 



(1) 



where 

+00 ^ 

W--=E(iT^ 

Under Vn , on T-L-r the process (^xj^ 
law to the Gaussian process (^xj^\ xj^\ 



„( 2 ) 



,(3) 



< +OO 



,( 4 ) 



N \ X^' N \ Xf^' N \ X^' N ^ converges in 



X^\xj^J solving the following linear 
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stochastic differential equations 
dx£\t) 



lh -h 2 )x£\t)~huX { *\t) 
lh -h 2 )xP(t) + hu J X^\t) 



dt - 
dt - 



T2 dW * ){t) 



dX<£\t) 



2 
2 



X%\t) + huX£\t) 



dt 
dt 



V2 



where S\h is Kronecker delta and 

• (x^(0),X^ ) (0),X^\Q),xj^ ) (0)j ^ is a centered Gaussian vector with Cov (x^ (0) , (0) 
for i =/= j or h =/= k and Var (x^(0) \ — | for any i, h. 

• wff\ Wff*^j are independent standard Brownian motions, that 
are independent of (x^ ] (0) , X^ ] (0) , xf ] (0) , X^ (0)) 

V / h>l 



Note that the randomness of the field appears only through the parameter ui in the 
dynamics of fluctuations. The only source of stochasticity is due to the Brownian 
motions. 

We now proceed to the analysis of the critical regime, i.e. for = 8 C A 2, where 
6 C is given in (l2~5j) and, under (HI), 8 C = 1 + 4u 2 . We make the following further 
assumption. 

(H2) w<%. 

Under assumptions (H1)-(H2), we have sufficient control of the spectrum of £, as 
operator in £ 2 ([0, 2n) x { — 1, 1}). In particular, £ can be diagonalized in the critical 
regime, as stated in next Lemma. 

Lemma 3.2. Under assumptions (H1)-(H2), ker(£) ^ if and only if 9 = l+4w 2 . 
In this last case the spectrum of £ is given by 



1 



0, -- + 2u/ } U 



Spec(£) 
with corresponding eigenspaces 



ker(£) 

Eig (-|+2o; 2 ) 

Eig ( + ikuj 

) I < 
Eig 1—4 — ikuj 



k 2 

— ±ikLo, k eZ\{-i,o,+i} 



/ (1) (2)' 

span \v\ ,v\ 

} (3) (4) N 

span \v\ ,v\ 

) (1) (2)' 

span lv k ,vl 

) (3) (4) N 
span \v k ,v k 



where 



Vi(x,r]) 
v[ 3 \x,n) 
Vk\x,n) 

vr (x, n) 



cos x — 2ujnsmx 
r\ cos x + 2uj sin x 
sin(fcir) — irjcos(kx) 
sin(fcir) + irjcos(kx) 



V^\x,T]) 
Vf\x,Tl) 
vr (X, 7]) 
V ( k 4) (X, 7]) 



sin x + 2ujn cos x 
2uj cos x — r\ sin x 
cos(fcx) + in sin(fcx) 
cos(fcx) — insh\{kx) 



(29) 
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In the critical regime 9 = 9 C = 1 + 4w 2 the variance of the processes 

U^ N Ht) := x[ hN \t) - 2coX[ 4 ' N) (t) and U^' N \t) := x[ 2 ' N \t) + 2coX[ 3 ' N) (t) , 

which are the fluctuations of the empirical averages corresponding to the directions 
generating the kernel of operator £, diverge as t — > +00. A sharper description of 
the large time fluctuations is obtained by considering more "moderate" fluctuations: 

p N := N~^p N . 

We will obtain asymptotics, as N — >• +00, for the signed measures pN(V~Nt). Note 
that these measures are completely characterized by their integrals 

vt N \t) := J v { ,:\x, V )dp N (VNt), (30) 

with h > 1 and i = 1,2,3,4. 

Theorem 3.2. Assume 9 C = 1 + 4w 2 , and oj < ^h^- Under Vn the processes 
(yJ? ,N) (tf) , for i = 1,2,3,4, and v} 3,N) (t),v} 4,N) (t) collapse to zero in the 

sense of Definition \2.1[ while the process (v^ N ^($),V^' N ^{t) \ converges weakly to 
the unique solution (V^'fy), V^ 2 \t)) of the stochastic differential equation 



dV^Ht) = - i il^J ) & ) V^(t) [(V^(t)) 2 + (V^(t)) 2 ] dt + J^dWV{t) 

dV^it) = -^S^^V^(t) [(V^{t)) 2 + {V^(t)) 2 ] dt+Jm^dW^tt) 

. V^(0) = V^(0) = 

where W^ 1 ' and are two independent standard Brownian motions. 

In the case ^= < u) < |, the process (V^ 1 - 1 (t), (£)) explodes in finite time; 

the convergence above holds for the localized processes: for every r > 0, the process 

(V^' N) (t A T N>r ), V} 2 ' N) (t A T Nir )) converges weakly to (t A T r ), (t A T r )) , 

where 

Tn,t ■= mf [t > : (v^ N \t)) 2 + (vf™®)' > r} 

T r := inf jf > : (V«(t)) 2 + (V (2) (i)) 2 > rj . 

By Theorem 13.21 we can derive the limiting dynamics of the critical fluctuations for 
the homogeneous model p = 5$. They can be obtained as a particular case setting 
u = 0. 

Theorem 3.3. Assume 9 C = 1. For h > 1 integer, let 

F fe (1,JV) (t) := /" cos{hx)dp N {VNt) and Y^' N) (t) := y shx(hx)dp N (VNt) . 

Under Vn the processes \ Y^ l ' N \t)^ , for i — 1,2, collapse to zero in the sense 
of Definition \KT\ while the process ^Yj 1,N \t),Y^ 2 ' N \t)j converges weakly to the 



16 



unique solution of the stochastic differential equation 

dYW(t) = -\YW{t) \(YW(t)) 2 + (Y^(t)) 2 ] dt+^dWW(t) 

dYW(t) = -{YW(t) [(yW(i)) 2 + (Y^{t)) 2 ] dt+^dWW(t) 

{ Y 1 (1 \0)=Y 1 (2 \0) = 
where W^ 1 ' and are two independent standard Brownian motions. 



4 Collapsing processes 

Before giving the details of the proofs of the results stated previously, we briefly 
present one of the key technical tool: a Lyapunov-like condition, that guarantees a 
rather strong form of convergence to zero of a sequence of stochastic processes. The 
first result (Proposition 14. 11) we state concerns semimartingales driven by Poisson 
processes, whose proof can be found in the Appendix of 5 . In the case where 
the driving noises are Brownian motions, the result takes a slightly simpler form 
f Proposition ^. 2 j) ; its proof is a simple adaptation of the one in [3], and it is omitted. 

Proposition 4.1. Let {£, n (t)} n >i be a sequence of positive semimartingales on a 
probability space (O, with 



dt„(t)=S n (t)dt+ / f n (t-,y)[A n (dt,dy)- A n (t,dy)dt}. 

Here, A„ is a Point Process of intensity A n (t,dy)dt on R + x W , where W is a 
measurable space, and S n (t) and f n (t) are s4 t -adapted processes, if we consider 
(^t)t>o o, filtration on s/, generated by A„. 

Let d > 1 and Ci constants independent of n and t. Suppose {n n } n >i, {a n }n>i 
and {/3 n } n >i, increasing sequences with 



Kit a, 



> 0, n n a n — 



>o, <%- +v 



> 



E 



(Cn(0)) d <C ia - d for all n. 



Furthermore, let {T n } n >i be stopping times such that for t G [0,t„] and n > 1, 
S n (t) < -K n 5£, n (t) + f3 n C 2 + C 3 with 5>0, 
sup \fn{t,y)\ < C A a~ Y , 



(f n (t,y)) 2 A n (t,dy) <C 5 



Then, for any e > 0, there exist C§ > and uq such that 



sup 2? \ sup £„(t) > C 6 (nn^ 1 V a n K n A 

n>n Lo<t<TAT,- " ' 



(al) 
(a2) 

(o3) 
(o4) 

(o5) 
(31) 



Proposition 4.2. Let {£n(£)}n>i be a sequence of positive semimartingales on a 
probability space (fl,s/, s ), with 

dUt) = S n {t)dt + f n (t,i)dWi{t) . 
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Here, (Wi)™^ are independent standard Brownian motions which generate a filtra- 
tion (&f t ) t >o> an d S n {t) and f n (t,i) are srf t -adapted processes. 

Let d > 1 and Ci constants independent of n and t. Suppose {K n }n>lj {a n }n>i 
and {fi n }n>i, increasing sequences with 



3 1 n-t+oo . n-H-oo „ -la n-y+oo 



Kn Ct„ 



J r\ — 1 

0, K n a n 



> 0, K n f3n 



±>0 



(61) 



and 



E 



(e«(0)) ] <C ia ~ d for all n. (62) 
Furthermore, let {T n } n >i be stopping times such that for t 6 [0,t„] and n > 1, 

S« (t) < -K n 6£ n (t) + /3 n C 2 + C 3 with 5 > 0, (63) 

£/«M 2 <C 5 . (W) 

i=l 

Then, for any e > 0, there exist C§ > and no suc/l t/iat 

sup ^ I sup £„(t) > C 6 (kIc^ 1 V a„K" lN ) 1 < e . (32) 

n>n iO<t<TAr n v 'J 

5 Proofs for the Random Curie- Weiss Model 
5.1 Preliminaries 

Proof of Lemma 12.11 An equilibrium probability density for (Q} must satisfy 



V 7 



= 0, 



which is equivalent to e P rT \ m *+ r i) (a, n) = e^ <T< -™ l * + ''- l (7*(— <7, rf), where m* is defined 
in ([5]). Solving, we obtain 

< 7 ,((7,r ? )= e ^+''>, 

with the normalizing constant 

Z,= f e^ {m ' + ^da = 2cosh(p{m,+r 1 )) , 

J. 5* 

and the proof is complete. 

Proof of Lemma l2.2l Obviously £ is a linear and continuous operator. We have 
to prove that, if 0i, 02 eL 2 (^),then $ s {24\{n))<h{y) v {d'n) = f a <h(v) i&faiv)) v(drf). 
Thus, 



(£0i fa)) MvMdn) 

cosh(/3(m» + rj))<t>i{rf) - P 
cosh(/3(m* + r)))(f)2(r/) - P 



s cosh(/3(rn* + 77)) 

4>2{rj) 
$ cosh(/3(TO* + 77)) 



(i(dr)) 



and the proof of self-adjointness is completed. 
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Proof of Lemma l2.3l To prove positivity of £ we have to show that if <f> £ l?(y), 
then fg (£<f)(r])) <j){rj)v(drj) > 0. Indeed we have 



(£^(77)) HvHdv) = 



cosh((3r])(j)(ri) — (3 



g, cosh(^?7) 



H{drj) 



<j>(ri)v{dri) 



— / cosh 2 (/3?7)(?!) 2 (77) £ p(dri) 

■ 1 KH " V cosh^y V 



cosh((3r])(j)(ri) f , K d v) ) > 0. 

cosh (prj) 



where we have used Jensen's inequality for the probability /3 ^j^g^y- Moreover, 
equality holds true if and only if cosh(/3r])(f)(r)) is constant; therefore the null space 
of the operator £ is generated by the functions of the form <j>{rj) = cosh 1 ^ ? ^ . 

5.2 Expansions of the Infinitesimal Generator 

Proof of Proposition I2TT1 By direct computation, and Taylor expansion of tp, 
we obtain 



L N ip ( / a<f>(r])dp N 

N 



.7=1 L \J 



N 



N 



3=1 



1 + V^7 (-@% ( °dp N \ +o(-^ T 



2<JA 



h=l 

2 



i=i 



,fe=i 



N 

n 

{ ^MVj^j cosh (P( m * + %■)) -sinh(^(m* 



N 



i=l 



3=1 



N 

{?_ 
2N 

P 3 
6N% 



N 



adp N J ^>i(»?i)[cosh(/3(m» + rjj)) - er,- sinh(/3(m» + rj ))] 

crdpjv I ^2 faiV^Wj cosh(/3(m* + - sinh(/3(m* + 7/j))] 
' 3=1 

J adp N ^j i (77 i )[cosh(/3(m st + rjj)) - crj sinh(/3(m* + 77^))] 



+ $W>0) { MVj)MVs)[<x*h(P(™* + Vj)) ~ °3 sinh(/3(m» + %))] \ + o(l) 

3=1 



i,fc=l 



We now represent all the terms as integrals with respect to the measure /5at; since 
-o-cosh(/3(m* +77)) + sinh(/3(m* + 77))] (j>i(r))q*(da, drj) = 



and 



[cosh(/3(m* + 77)) — <r sinh(/3(?7i* + 77))] cf> i (r])q ill (da,dr]) = I 



0(»7) 



^ cosh(/3(m» + 77)) 



H{drj) . 
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we obtain 



L N ip / a(j}{rj)dp N 



cosh(/3(m* + rj))<j>i(rj) - ft 



@ cosh(/3(m* + T])) 



dp 



N 



sinh(/3(m* + rj))4>i(i])dp N 



N 
f?_t 
2N \ 

6N 
J*_ 

6iV§ 



adpN / [cosh(/3(m* + rj)) - crsinh(^(m* + rj))](j)i{ri)dpN 



adpN ) J [o cosh(/3(TO* + 77)) — sinh(/?(m* + r]))]<fri(r))dpN 
3 



erd/5 



@ cosh(^(m* + rf)) 



p,{drj) 



adp N / [cosh(/3(m* + 77)) - a sinh(/?(m* + T]))]<j>i{r))dpN 



2 £ 9^(0 {/ 



cosh(/3(m* + 77))' 



M(*7) ^ + o(l) , 



from which (fTU|) follows. The fact that the remainder o(l) has the form (TT2"j) and sat- 
isfies (|13[) is implied by the Lagrange form of the remainder of the Taylor expansions 
we have used. 



Proof of Proposition 12.31 It is obtained by a simple rescaling of the last ex- 
pansion of L^-tp (J a<p(rf)dpN) seen in the proof of Proposition 12. II The details are 
omitted. 



5.3 Collapsing Terms 

For N > 1, M > define the family of stopping times 



inf 

t>o 



Y} N \t) 



> M for at least a value of i = 0, . . . , m — 1 >, 



where the 5^ s have been defined in (TT5|) . In the rest of this section, we often 
consider the time-rescaled infinitesimal generator J/v = N^L^, where Ljy is given 
by (TH|) . Whenever we write 



we mean 



JjvV' / <7ip dpN, / aifxdpN, 



<J(p m -idp~N 



p N =p N (Nit) 



We later consider, for j 6 y and k & 2>, the counting process A^(j,k,t) which 
counts the number of spin flips of spins at such that at = j and rji = k, up to time 
N*t. We consider the following semi- martingale decomposition 



'{Y} N \t)) 2 = J N \( 



Y : 



(N) 



(t) dt + dM% Y 2 , 



(33) 
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with M.* N Y 2 the local martingale given by 



where we have defined 



(Y! N \t)j 



Yl N \t)-j 



2<pi(k) 



Ni 



A a N (j,k,ds) : 



(34) 



(35) 



and 



A%(j,k,dt) :=A%{j,k,dt)-m 



A(j,k,mt) 



-Pj(N 4 / adp N (Nlt)+k) ^ ^ 



■—\"{j,k,t) dt 

The quantity A^(j,k 7 dt) is the difference between the point process A^(j,k,dt), 
defined on y x f x l + , and its intensity X a (j, k, t) dt. The quantity A(j, k, Nit) 

indicates the number of sites % that at time N*t have a.- L — j and r/i — k and it is 
given by 



A(j,k,mt) 



N 

T 



i + 



i 



kNi 



rjdp N {t) 



J 
Ni 



adp N (t) 



| f^r y <?vdpN(t) - J^i]ta,iih(/3ri)p(dr]) 



(37) 



Remark 5.1. If we call (.At)t>o the filtration generated by A^, then the processes 



N 



and V 



CO 



are ^—adapted processes. 



For every index i = 1, . . . , m— 1, the following result holds. Note that it is stronger 
than the collapse of the processes \ J , in the sense of Definition 12.11 



Lemma 5.1. Fix d > 2, and assume the assumptions of Theorem ] 2. 3\ are satisfied. 
Then, for every e > there exist Nq such that for every M > there is a constant 
Cq > for which 



sup P < sup 

N>N [o<t<T/\: 



(yW(t)) 2 >c 6 iv-Hi-§) 



< e. 



(38) 



Proof. The main tool is Proposition ^. II However, some assumptions in Proposition 
14.11 are not satisfied uniformly in the environment. We therefore will condition on 
the event 



A K ■= { v e 



smh(J3r])ipi{r])dp N 



cosh(/3r])ipi(r])dp N 



< K 



The random field r\ is i.i.d., so it satisfies a standard Central Limit Theorem. There- 
fore, we can choose K > such that for every N > 1, 

*W < §• 

Constants below are allowed to depend on K; this dependence is omitted. We are 
left to show that, for every M > there is Cq > such that 



sup P K 

N>N 



sup 

0<t<TAr" 



( Y ( N \t)) 2 >c 6 N-H^) 



< 



(39) 
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where Pk{ 1 ) : = P{ • |AkO- To prove we check the conditions in Proposition 
EQ 

Step 1. We set kn :— N*, ajv ■— N» , (3 N = 1. Clearly (jaTj) in Proposition 14.11 
holds. 

5iep 2. We check (|a5) of Proposition [O] i.e. 



(*f°(o) 



2r?" 



<dN~^ for all TV. 



(40) 



We start noticing that a Central Limit Theorem applies to the processes J aipi(rf) dp^ifi), 
since the random variables (ctj(0), V'i( 7 7j))j=i are independent; so, in the limit 

as N —> +oo, N^Y^ N \o) converges to a Gaussian random variable and, since 
(oj(0)Vt(?jj))|Li are bounded random variables, there is convergence of all the mo- 
ments. In particular (|4U| holds. 
Step 3. We check Jo3j of Proposition ^. 11 i.e. 



Jn 



Y 



(AT) 



(t)<-7V^(y.W(t) 



C 2 



(41) 



for suitable constants (5, C2 > 0, which are allowed to depend on M, and all 



t€ [0,Ttf] (we recall that j3 N = 1). Letting X 



._ £ 



J adpN, we write 



cxp[±X] = l±X + R ± . 



Using this expansion we can perform the computation as in Proposition 12.31 but 
keeping track of the remainders : 



Jn 



(yt N) y 



N 



N 4 [cosh(/3?7j) — tjj sinh(/3?7j)] exp 



Y- 



(N) 



V [cosh^) - a j sinh(^)] ( 1 - / trdfrv + 

r - , \ N* j 



-47V3F- 



(AT) 



-4Y> 
smh(/3r))ipi(r])dp N 



Ni 



a&tfidpN + AY> 
-4r/ JV) /3 / adp N / a smh(Pr])(pi(rj)dp N 



+ 4jrW__ f ,jdp N [ cosh(pr])<pi(r])dp N 
Ni J J 



N 



T7T H [ cosh 0%) + a j sinh(/3%)] p?^) 



iV3 



J'=l 



4£ 
ATI 



AT 



crdpAr ^2 [cosh(/3^ ) + oj- sinh(/3?7 i )] ajipj (^) 



AT 



V] [cosh(^?7j) + ctj sinh(/3%)] 
j'=i 



-4K 



iV5 



(42) 
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The first term of this last expression is 



We are left to show that all remaining terms are bounded, for t 6 [0, rj!f], rj € Ak 
and assuming that in (j42]) . is evaluated at time Nit. We immediately have 



Y} N) (t) <M, 



smh(Pri)ipi(r])dp N 



cosh(/3r])ipi(r))dp N 



< K. 



All remaining terms in (1421) are of the form 

af( V )dp N (Nh), 



for some real valued /. Since (fh^—o form a basis for the vector space of these 
functions, we can write 

m— 1 

/ = E ahLph - 

h=0 



Thus 



af( V )dp N (N*t) 



m— 1 



< 



E ki ^ w (*) 



/i=0 



< CM, 



where C depend on m, on the combinators an, but not on N. As a consequence 



|i?±| <sup<je* : |*| < -^r 

N •* 



With all this, ((42]) implies 



A? 



adpN 



< -4A 



2AN 



cr(i/5 



A 



< 



f3 2 M 2 
2N^ 



m (y} n) Y + c(M) 



for some M-dependent constant C(M). 

Step 4- We check (fa4")l of Proposition l4.il i.e. (see equation (IMIO 



sup 

u>en,j£.y,t<T* 



7 U) 



yW(t) 



< c 4 iv- 



(43) 



For i < rjy , we easily have 



V (j) 



Ni 



Ni 



< -5-(l + M) sup{^(fc) + |^(fc)|} < C 4 JV-s 
A* fee® 



5iep 5. We check ()a5|) of Proposition ^. 11 i.e. (see equation ([36 



E 



A CT (j,fc,t)<C 5 



(44) 
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— (j) 

Recalling the definitions of V 
(|3"S")) and in (f3"6")> . we have 



E 



V 0) 



A CT (j,M) 



and X(j, k, t), which can be found in 



A* \Mj,k,Nh)\ 



2<Pi(k) 

at! 



Boundedness of this last expression for t G [0, rjjf ], r\ € Ajc follows readily by bound- 

ednessof y} N) (t) and / adp N (N*t) (see step 3), and the fact that \A(j,k,Nh)\ < N. 
Step 6. Conclusion. It is now enough to use (j3i~j) . □ 



The next step is to prove, for every e > and N > 1, the existence of a constant 
M > such that 

P{r^<T}<e. (45) 

This fact, together with Lemma l5.lt implies the processes Y{; N \t), . . . 
converge to zero in probability, as A grows to infinity, for t in the whole time 
interval [0,T]. As in (|3T))) . we can replace P by Pk for a sufficiently large K. 

The idea is to consider a martingale decomposition as in (|33p for V> ^o"') > wnere 

?A G C 1 has bounded first derivative, and is such that |af| > M implies ip(x) > M; 
for instance, ip(x) — \f\ + x 2 . We obtain 

1> (*f°(t)) = i> (Y (N \Q)) + f J N ^ (V W ) {s)ds + M t N ^, (46) 



where 



E 



A * 



-^(r w ( s ) 



A a N (j,k,d8) (47) 



with as in (f3"6"l) . The point now is to get bounds on Jn"4> ^q"')- We proceed 
as in (|42[) : the only difference is in the "gradient term", which is now 



^ (V w w - ^^P) - * or^) - (>r (*)) - 

with TZn < Proceeding as in (|4"2")l . it is easily seen that 



A3 



JniP (V W ) = 2ty' (V (JV) ) / wahiprfitpoWdpN 

- 2iP' (y (JV) ) p j adp N [ asinh^MvWN + O m (a*) , (48) 



where Om f-^') includes all term that, for t < rj^, are bounded by 
m 

20' (V (A,) ) / sinh(/37 7 )^ (7 7 )dp Ar 



absolute value of the term 



^1. The 

JV3 
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is bounded by CK , since ip' is bounded, and r\ € Ak- For the term 

one should notice that sinh(/3rf)(po(ri) is orthogonal to <fio(r]) in L 2 (v). This implies 
that / tr sinh(/377)(/Jo(?7)(i/5Ar is a linear combination of Y±\ Y^ , . . . , Y^\. Due 
to (j3"5j) . we can choose a constant C{M) for which this term is bounded by 

C{M)N~-^ 1 -i) 

for t < Tjm) with probability greater that 1 — | . Denote by B £ the event that this 
bound holds true. Putting all together, we have therefore proved that in Ak H B e 
and i < Tjfjf, 



C(M) 



This means that, by (|46[) . the inequality 



sup 

Q<t<T/\rM 



> M 



implies, for N and M large enough, that either 



y W(o) 



> cM 



sup I A4* 



for some c > 0. 
Thus, 

K' < r} c 



sup 

0<*<TAt? 



0<t<TAr? 



> cM 



Y 



(N) 



\(t)\} 



> M 



m—1 

1=1 



sup 

0<t<T/\r* 



Yf N \t) > m| u { r (Ar) (o) > cm} u 



u 



sup | |>cM^U4^UB e c 

0<KTAt" 



and we obtain the following inequality for the probability of the interested set 

m — 1 



P{t™ <T}<-e+Y, P \ SU P 

i=l •> 0<t<TAr£ f 

•{|y W(o)| > cM } 



> M 



P 



sup |M^| > cM J> . 

0<*<TAri;/ 



We estimate the three terms of the right-hand side of the inequality. 
• For every i = 1 , . . . , m — 1, thanks to (1381) we have 

e 



P 



sup 

0<t<TATj? 



Y{ N \t) 



> M > < 



12 



where for M large enough. 
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• Since at time t = the spins are distributed according to a product measure, 
Yq N \o) is times the sample average of independent, bounded random 
variables of mean zero. Therefore, for some constant C > 0, 



E 



y o w ( o) 



< 



c 



and in the limit as N — > +oo, we have convergence to zero in L 1 and then in 
probability. Therefore 



P{|F W (0)|> cM} <± 



for N sufficiently large. 



We reduce to deal with E 



{^4 



in fact, by Doob's maximal inequality 



for martingales (we refer to Chapter VII, Section 3 of [T2]) we have 



E 



sup \M* N J >cM}< 



0<KTAt" 



It is therefore enough to show that E 



{cMf 



is bounded uniformly on 



TV and M. By (|4"T)) and since tp is Lipschitz, we have (see also (|3"6")l ) 

C 



E 



(M 



< —^E 



jeJ",kes> 



(49) 



Since, by (13^1) . X a (j,k,t) < CN* for some constant C, the boundedness of 



E 



is established, and the proof of (f4"5j) is completed. 



5.4 Identification of the Limiting Generator and Convergence 

We are going to show that, in the limit of infinite volume and t € [0, T], the process 
Yq (t) admits a limit in distribution, that we will be able to identify. 

First, we need to prove the tightness of the sequence |lo^(*)} • This property 

implies the existence of convergent subsequences. Secondly, we will verify that 
all the convergent subsequences have the same limit and hence also the sequence 

I v(")w I must converge to that limit. 
L J N>1 

Lemma 5.2. The sequence ^YQ N \t)} is tight. 

Proof. Following [5], we use the following tightness criterion: 

a sequence of processes {£Ar(i)}jv>i on T>[0,T] is tight if 
1. for every e > there exists M > such that 



N 



SUpPl SUp |£y(t)| >M \ <£, 



te[o,T] 



(50) 



2G 



2. for every e > and a > there exists 8 > such that 



sup sup -P{|6v0"2) - £jv(ti)| > a} < e . 

N 0<t 1 <t 2 <(t 1 +6)AT 



(51) 



where the second sup is over stopping times t\ and t-i , adapted to the filtration 
generated by the process £jv- 

We must verify the conditions ([5U|l and ([5"Tj) hold. Since we have already proved 
that for every e > the inequality P{t^ < T} < e is true for M sufficiently large 
and uniformly in N, it is enough to show tightness for the stopped processes 

We have already shown that, for M large enough 



P ■ 



sup 

0<t<TAri 



y W(t) 



which yields (j5"0)l . To obtain (|5"Tj) . we notice that 



(52) 



where we have denoted 



J l N,Yo 



2 r T 2 „ 



and A^r is as in (|3T>1) . As in the proof of Lemma [5TT1 one shows that both Jn (Yq N ^ 
and the quadratic variation of M t ^'y are uniformly bounded in N, from which (|51[) 
follows for the processes ^Y$ N \t A t^ 7 ) j . □ 

Lemma l5~2l implies that there exist convergent subsequences for the sequence |Yo^(*) j 
With abuse of notation, let I Yq (t) \ denote one of such a subsequence and let 

^ J n>l 



ip £ C£. The following decomposition holds 



^(F W(*))-^(y o Cn) (0))= j J^{Y^\u))du + M 



L n,'ip i 



(53) 



where 



J n ^{Y^ n \t)) =21/ (V ( "V)) Ini j t a nh(j3 V )dp n (t) 

+ p / ad~ Pn (t) / dp n {t)-/3 / adp n {t) / atanh(pr))dp n (t) \+o M {l)- 



The remainder Ojvf(l) goes to zero as n — > +oo, for t < r„ . If we compute the 
limit as n — > +oo, using the facts that a Central Limit Theorem applies to the term 
J taiih(j3ri)dp n (t), the integral J dp n (t) is zero since p n is a centered measure, and 
the process J a ta,nh((3r])dp n (t) collapse since tanh(/3?7) and ^Po(t]) = cosh 1 ^ ]? ^ are 
orthogonal in L 2 (y), we have, in the sense of weak convergence of processes: 



(Y Q {n) (t A rf )) ^±22, j^y^ A ^ with j^y^t)) = 2 Jf (y (t)) 
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and where Jrf? is a Gaussian random variable. Then, because of (f5U)) and (|4"5)) . we 
obtain 

Mi ti> M% := ^(y (t)) - V(ib(Q)) - r JHYo(u))du , 

™ Jo 

for t G [0,T]. We must prove the following Lemma: 

Lemma 5.3. Jvl 1 ^ is a martingale (with respect to t); in other words, for all s, t G 
[0,T], s < t and for all measurable and bounded functions g(Yo([0, s])) the following 
identity holds: 

E[M%g{Y {[Q, s]))] = E[M%g{Y Q {[Q, «]))] . (54) 

Proof. We begin by showing that ((54)) follows from the fact, that will be proved 
later, that for every t fixed, {M.^ ^}n>i is a uniformly integrable sequence of random 
variables. 

Since M. n ^ is a martingale (with respect to t) for every n, we have that for all 
s,t £ [0, Tj, s < t and for all measurable and bounded functions <?(>o([0, s])) 

£[M^ 5 (Yo([0,s]))] = S[M, s v0 .g(r o ([O, S ]))]. (55) 

Now, as we have seen, .M^ ^ and A^^ ^ have a weak limit; this, together with 
uniform intcgrability, imply convergence in L . Thus (|54p follows by taking limit in 
(ESJ). 

It remains to check that {M.^ ^] n >i is a uniformly integrable family. A sufficient 
condition for uniform integrability is that sup„ -E^A^ ^| 2 ] < +oo (see again [15]). 
This, however, is exactly what we have done already in (|4T)|) . □ 

Proof of Theorem l2.3l We have shown that any weak limit of Yq ( • ) solves the 
martingale problem with infinitesimal generator J, which admits a unique solution. 
It follows that all convergent subsequences have the same limit and so the sequence 
itself converges to that limit. 



6 Proofs for the Random Kuramoto Model 

Throughout this section we assume u> < ^tjj even though this assumption will be 
relevant only starting from Section [6.31 Whenever needed, we will comment on the 
necessary changes to cover the case ^= < w < \. 

6.1 Preliminaries 

Proof of Lemma 13.21 If ip(-, •) belongs to the null space of £, then Zip = 0. 
Therefore, we require that 



We solve the ordinary differential equation (|56|) . Having defined 

A := ^- J cosy<p(y,<;)q*(dy, efc) and B := ^- J sin?/^(y,<r) q*(dy,d<>) , (57) 

the solution is ip(x, rf) — 2{B — 2Auji]) sin a; + 2(A + 2Buji]) cos x; this function yields 
a solution of (|5rJj) provided that it satisfies the self-consistency relations ([5T|) . but it 
does for every value of A and B. Then the two directions which generate the kernel 
are sin a: + 2uir]COSx and cosx — 2uir]sm.x. 
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Remark 6.1. In the case that 9^1 + 4w 2 , the unique value for which the self- 
consistency relations in (l57l) are satisfied is A = B = 0, meaning that at the critical 
point the kernel of the operator £ is two-dimensional, while it is trivial for all the 
other values of the parameter 9. 

The part of the statement of Lemma 13.21 concerning spectrum and eigenspaces is 
easily proved by direct computation, and the fact that the set {vu : k > l,i = 
1, 2, 3, 4} spans a dense subset of £ 2 ([0, 27r) x { — 1, 1}). 



6.2 Perturbation Theory 

In the rest of the section, we often consider the time-res caled infinitesimal generator 
Jn = \/NLn, where is given by f|27|) . To determine the limiting generator J, 
we need to apply the first order perturbation theory. The methodology for treating 
a perturbation problem has been developed in the paper |17) and extends the earlier 
works done in []3J [TB] . 

It will be useful to keep in mind the following simple fact, which is just a restatement 
of Proposition 13. II 



Proposition 6.1. Under the assumptions of Proposition^^ we have 

Jn^ {^J 4>{x, v)dp N ^j = VNL (1) 4> + N^L {2) ^ + L (3 V + N-iL^ip (58) 
where 

L {1) ij) := djip (J(j>(x,r])dp N J J £<f>i(x,r))dp N 

L (2) i/> := e^djip yj (j){x,rj)dpN^ J '^^'^ sin(y ~ x)dp N dp N 

L (3ty := 1 £ ^ y <p(x, V )dp N ^ J ^Lfan^fatfdq. 

L[A) ^ := \ ]C d ik^ (^J <P(x,v)dpN^J J ^-{x,T])^-(x,T])dp N . 

As first step fSection UP)) we show that for every cf> e span jwj[ 3 \ v^\v^ : k > 2, i = 1, 2, 3, 4 j, 
the process 

{x,r/)dp N (VNt) 



collapses to zero in the sense of Definition 12.11 We are therefore left to understand 
the behavior as N — » +oo of the two-dimensional process 



(v^ N \t),v} 2 ' N \t)) ■■= [J vPd?s(VNt),J v^dp N (VNt) 
For this reason, for ip £ C 2 (R 2 ,K), we need to control 

Jn^ (^J dp N , J v[ 2) dp N 

The first term in the r.h.s. of (15"5|) vanishes, since v± , vf^ £ ker(£). In order 
to compensate for the second diverging term N^L^tfj, one introduces a "small" 
perturbation of ip of the form 

ip N = ip + N~iipi, (59) 
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for some ipi to be chosen. We obtain 



J N ip N = N* 



L (i) i> + L (1) i/>ij + L (3) V + i (2) V>i + o(l). 



(60) 



In order to avoid divergence, tp\ should be chosen in such a way that L^ip+L^tpi = 
0. At a purely formal level we are led to set 



in :=-(l«)~Vv, 



which gives 



JnVn > 



-.(3) _ r.(2) ( L W 



L 



(2) 



ip =: Jip. 



(61) 



(62) 



The operator J is therefore the candidate for the generator of the limiting process 
(yW,y( 2 )). In order to make a rigorous proof out of this formal argument, the 
following two steps are needed: 

1. The operator (L^) 1 has to be properly defined. 

2. From the above convergence of operators one must derive weak convergence 
of processes. 

Step 2 will be dealt with in Section |6~^1 through standard martingale techniques. 
We consider now step 1. The needed computations are rather long, but follow few 
basic ideas, that we now illustrate. First observe that 



L {2 H(J v^dpvJ v{ 2) dp N ) =8j2dii>(-,-) f 



dx 



■{x, -q) sin(y — x)dp~NdpN . 

(63) 

We give the details for the term J v) sin(y — x)dp~Ndp~N, the other being 

similar. Letting 

y^ N] ■■= J v?dps, 

by applying standard trigonometric formulas we obtain 

a (!) 



dx 



■{x, rj) shx(y — x)dpNdpN 



4(1 - 4w 2 ) 



iuj 



2(1 - 4w 2 ) 



(4,jV) _ y(2,N) 



V^ N) -2uV^ N) ) (V 2 (1 < 



N) (3,N) 



(64) 



This means that L^^fjv^dpNj v { dp~N I i 
the form 



is a linear combination of terms of 



ditp (^J dp N , J v{ 2) dp N ^j J v[ 3) dp N J v { 2 h) dp N , (65) 
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i = 1,2, j, h = 1,2,3,4. If we denote by X 3 k the eigenvalue of £ corresponding to 
the eigenfunction \ in the critical case 9 = 1 + uj 2 , we easily obtain 



—j ~~^ di^> (^J v^dpN, J v[ 2) dp N ^J J v[ 3) dp N J v 



2 dPN 



dii> (Jv^dp N ,Jv^dp N J J v[ 3) dp N j v 2 h) dp N ; 



this defines 1 for the whole expression in (j64|) . Thus, the perturbation (l6Tj) 

is now well defined. 

A further comment is relevant. In the expression for the limiting generator in (|62p . 
the quantity 

appears. Moreover, we have seen that (L^) 1 L^ip is linear combination of terms 
as in (p5|) . We will prove later that, when evaluated at time y/Nt, 

• the sequences of processes J v^'dpff, j = 3,4, and J v^" dpNi h = 1,2,3,4 
collapse to zero; 

• the sequences of processes J %yf'dpn, j — 1,2 are tight. 

In particular, the processes (L^) 1 lP^i\> collapse to zero. We then have to apply 

L (2) 

again. It is easy to show what follows. 

• When dit/j (^f dp f dp n^J f dp n f v^dpN has j = 3, 4, i.e. it has 
"two collapsing factors" , then 



L 



(2) 



diip ( / v^dpN, I v^'dpN I / v\ J 'dp N I v^'dp 



,(2 



is still collapsing to zero. 
• When j = 1,2, non collapsing terms in the expression above, arise from 

diip (J v[^dp N> J v^dpN^J J v[ 3) dp N J (x, n) sin(y - x)dp N dp N , 

since when the Prostapheresis formulas are applied to g x (x, 7?) sin(y — a;), 
terms of the form J v^dpjy J dpN, j, I = 1,2, appear. 

Carefully performing a long but straightforward calculation, one obtains the follow- 
ing statement. 

Proposition 6.2. Up to collapsing terms (as the symbol ~ is intended to mean) 
we have 



-.(2) 



(l + 4^ 2 ) 2 (l-8^ 2 ) TA(1;jV) 
'4(l-4w 2 ) 3 (l+w 2 ) 

(l+4^ 2 ) 2 (l-8^) T ^ (2|jY) 
4(l-4w 2 ) 3 (l+w 2 ) 



with 



V, 



/ (i) 
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6.3 Collapsing Processes 

From now on we always assume = 1 +4cj 2 , with lu < i. In what follows, it is more 
convenient to work with the following real- valued basis of L 2 ([0, 2ir) x {—1,1}): 

{v?\yjp :i = 1,2,3,4, h > 2}, 

where 

y^P(x, rj) := coshx y^ (x, 77) := sin hx y^ (x, 77) := r) cos hx y£\x, rf) := rjsmhx. 
We also set 



Y, 



and write Y<f' N) (t) for J yjpdpft(VNt). 
For r > 1 define 



y h dpN, 



WPN 



(F <=,») ) % (F <«,», ) %^ s _j_ (r <. 



N) 



(66) 



Clearly, showing that the sequences of processes J v±dpN, j — 3,4, and J v^dpN, 
h = 1,2,3,4 collapse to zero is equivalent to show that the sequences of processes 
J v± dpN, j = 3,4, and J y^dpx, h = 1,2,3,4 collapse to zero which, in turn, is 
implied by the fact that the sequence ||/5jv||^ collapses to zero. All processes here 
are meant to be evaluated at time y/~Nt. For N > 1, M > define 



r*f := inf 



p N (VNt) 



> M or 



V^ N \t) 



> M or 



V^ N \t) 



> M 



Our first result concerns collapsing of the stopped process 



Lemma 6.1. Fix d > 2 and r > |. Then, for every e > and M > 0, there exist 
No > and C5 > 0, for which 



sup P sup \\p N (VNt)\\l > C 5 N^-^ \ <e. 



N>N„ 



0<«TAt" 



(67) 



Proof. We apply Proposition ^. 21 We set kjv = V^V, ctN = N 4 , @n = N 4 . Condi- 
tions (fFTj) and (|b2|) of Proposition l4~2l are easy to check. We are therefore left to check 

2 

conditions (fb3j) and (|M|) . We observe that p^{\/Nt) admits the semimartingale 

r 

representation 

_ 2 _ N d 

p N (VNt)j r = J N \\p N f r (VNt)dt + NiY,g^ \\pn\\1 dW.it), 

3=1 J 

where {Wj(i) : t > 0,j = 1, . . . , N} is a system of independent standard Brownian 
motions on [0, 2tt]. We show the following inequalities for every t <E [0, rff], which 
imply (O and ([55]): 

Jiv||/5iv||^(VFt) < - Q- 2 ^ 2 ) VA 7 ||/5 J v||'(\ / A 7 t) + C , A^3, 



JV 

3=1 



^ „- „2 
PJVIIr 



< C 



(69) 
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for some constant C, that is allowed to depend on M. 
Step 1: proof of (|68]l. We use (f58|): 



^iv||Piv||' = V^i^Ml^ll' + ^^^Ml^ll' + ^Ml^ll' + ^^^MlPivIl'. (70) 

We begin to deal with \\f>N\\r- Due to uniform convergence of the series defining 
ll/5?v|lr' we can apply term by term. For i = 3, 4 

(| i^dp* ) 2 = - (1 - 4^ 2 ) (| w«dp* ) 2 . (71) 



Also, by direct computation, 

4 



£ (1) E(/^) = -^ 2 E(7^ 

i— 1 ^ ' i— 1 ^ 



Letting A := 1 - 4w 2 > 0, by {7T]) and {72]) we obtain 



i (1) ll/5i V ||" = -A||piv||" + 



/t>2 v ' i=l v,y 



(72) 



(73) 



We now compute L( 2 )||pjv||J!- A "typical" summand with h > 2 of the infinite sum 
giving \\p N \\l, is 



L 



(2) 



y^d^l = 2(1 +4^) / yfdp 



>N 



dy (l) 

— sin(y - x)dp N dp N . 
ox 



By using Prostapheresis formulas, one realizes that J y^dp~N J sin{y—x)dpisrdpi\r 
is a linear combination, with uniformly bounded coefficients, of terms of the form 



v (j,N) v (i,N) v (i,N) v (j,N) 



(74) 



bumming over ft > 2 and observing that, for t G [0,r^], Y^' N \t) < cM for some 
constant c, we obtain (omitting the evaluation at y/Nt) 



(75) 



for some M-dependent constant C(M). As far as the first two summands of ||pjy| 



are concerned by similar arguments a rough bound for t G [0, tM\ of the form 

2 



E ( / W*p» 

i=3A \ J 



is obtained. Putting together ([73|) . (ff5|) and ([76 

ViV^^llpjvll' + TVU^IIpjvll 2 

< -VNX\\p N \ 



< C{M) 



(76) 



h>2 y ' i=l yJ 

Nic(M) J2 tj-^w E (/ y [ h d p 

h>2 y ' i=l w 



< — v NX WpnWi- + N^C(M) 



(77) 
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where the last inequality holds for N sufficiently large so that 

VN (A — 2h 2 ) + C(M) < 

for every h > 2. 

Consider now the term L^> \\pN\\ 2 r - We have 



dv{ 



(4)' 



dx 



dq* 



2 J2 r 1 + h 2yi,J 

h>2 v ' i=l J 



dx 



dq*. (78) 



By the simple bound 



using the fact that for r > 



dx 



h 2 



h>2 



< +oo, 



from (j78|) we get 

£ (3) ||piv||J!<C (79) 

for some constant C. The treatment of the term L*- 4 -* ||/3at||^ is quite similar, since it 
is obtained from (|78[) replacing q. M with p^v- Having pjy total variation ATi, we get 



L^\\~PN\\l<cm. 

By (EZl), dZll) and (EOl), (E3) follows. 

Step 2: proof of (169p . Consider the summand 



V ( h d PN 



of ||p»jv|lr- The summands containing are dealt with similarly. We have 



(80) 



d 
dxj 



Vh dpN 



2 

Ni 



N 



<9l 



■ (Xj , T}j ) , 



so that 



Vh dPN 



< 



2h 



Vh dPN 



Thus 



J h>2 K 1 i=l W 



< 



2 ^ — > /i ^ — > 



< 



2C 



7i>2 V 7 i=l X 7 
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with C 2 := ^2 h (i± h 2y t where we have used the Cauchy-Schwartz inequality. Sum- 



ming up, for t < t m 



M 
N i 



, , C 2 , C 2 M 



d Xj llrJV "V " N i ^"»r- N s , 

from which (|6l?f follows. □ 

Remark 6.2. For later use, we observe that the M-dependence of the constant C in 
(|6"5)l comes form the estimates in (|74p and (jTS")) . where the factor yP' N ' is estimated 
by a constant C(M). If we replace such estimate with the trivial one 



U,N) 



< N*. 



we obtain the following estimate, which does not require any stopping argument: 



J N \\p N \\*(VNt)<-(±-2^-} v.v:>xi!:'!v XI) -:- C\ X. 



which implies 



sup E 

JV>l,t>0 



\pn\\1 



< +00. (81) 



We can now prove the main result of this section, corresponding to the first part of 
Theorem [321 



Proposition 6.3. For every T > 0, the sequence (\\pn\\ t {^/Xt)j collapses 
to zero. 

Proof. Given the result of Lemma 16.11 all we have to show is that for every e > 
there exist M, Nq > such that 

sup P{t^ >T}<e. (82) 

N>N 

Consider the function 



%p(x,y) := yl + x 2 + y 2 . 

Note that tp has uniformly bounded partial derivatives, and ip(x,y) > min(|x|, \y\). 
We begin by observing that 



W f >T}c sup \\p N \\ 2 r (Vxt)>^-\ 

[0<KTAt*' z J 



u{ sup ^(^(t),^(i))>fi. (83) 

[0<t<TAT™ V J 

By Lemma P6 . 1 1 for AT large the probability 

sup \\pN\\l(VNt)>%) 

[0<KTAt* z J 

can be made arbitrarily small. Thus, ([52")) follows if we show that for every e > 
there exist M, No > such that 

sup pj sup ^(v^ N \t),V^ N \t))>^)<e. (84) 

N>N 0<«rAr« V ' 1 \ 
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For the proof of we consider the perturbation 

ipN = tp + N~4lpi 
as illustrated in Section HOI with 

As seen in Section 16.21 is a linear combination of terms of the form 

dii> (Jv i [ ) dp N ,Jvf ) dp N ^ J v[ j) dp N J v^dp N , 

and therefore, up to time Tpf , can be bounded in absolute value by some M- 
dependent constant C(M). This implies that, for every given M and for large 
enough N 



P 



sup ^(v^ N \t),V^ N \t))>^)<p{ sup ^(-)> 

:t<TAr« V J {0<t<TAT$f 



M 



By abuse of notation, we write ipN(t) in place of 

1>* (v r h (<,A °(t):/i = l,2;i=l,2,3,4) 
Consider the semimartingale representation 



(85) 



1>*{t) = ^jv(O) + / J N ip N (s)ds + M N (t), (86) 
Jo 



where 

M w (t) =N*pJ*^(s)dW j (s) (87) 

where {Wj(i) : t > 0,j = 1,... ,N} is a system of independent standard Brownian 
motions on [0,27r]. We have 

P\ sup ^j N (.)> K<p\^ N (0)> 

[0<KrAr« 6 J I 9 

+ P{ sup Jjv^Ar(t) > ^ \ + P \ sup M N (t)>^-\. 

[0<t<T/\T™ 9J J ^0<t<TAr^ y J 

The term P{i/V(0) > 4^} is easy to control, since the random variables (0) 
converge to zero in probability. We are therefore left to show that the probabilities 

f Ml 

P i Sup J N tpN(t) > — 



and 



Pi sup Xjv(0>^i (89) 



0<KTAt?; 1 



9 



are small for N large enough. We begin to deal with (J5SJ). By (|5U)l and the choice 
of tpi, we have 

J^Vjv =i (3 V + i (2 Vi+o(l), 
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where the term o(l) is bounded by C j^P for some a > 0. Moreover, it is easily 
shown that L^tp is bounded uniformly in N and M. To deal with L^ipi, we use 
Proposition 16. 2[ which gives 



= (l + 4c 2 ) 2 (l-8c 2 ) AUN) ( .s 
4(l-4w 2 ) 3 (l+w 2 ) 1 U 
(l + 4u, 2 ) 2 (l-8u; 2 ) T , (2 ,jv) 



4(1-4w 2 ) 3 (1+cj 2 ) 



-yr '(t) 



v! hN \t))\(v^ N \t)) 



W(v) 



collapsing terms, (90) 



where, again, the "collapsing terms" are bounded by ^rp- Observing that 



and since, by assumption, 1— 8a; 2 > 0, the non-collapsing part of (|9"0")) is nonnegative. 



We therefore conclude that, for t < rfj 



Jn*Pn < C + 



C(M) 



with C independent of M, N. This implies that the probability in (|88|) is arbitrarily 
small for M (first) and N (then) sufficiently large. 
We now deal with ([59"]) . By Doob's Maximal Inequality 



P< sup M N (t)>^-}< 

0<*<TAt" 9 



E 



(M/9) 2 

(A//9) 2 

Up to term bounded by C j^P , we can replace with -S^ in (|9ip . Moreover 



(91) 



dtp 



i=l,2 



1 9vi 
ATI <9x 



< 



C 



(92) 



for a constant C independent of M,N. Inserting this in (|9ip . we have, for some 
C, C(M) > 0, 



Pi sup XatW>^1< C 

0<t<TAr" « 



C(M) 



(M/9) 2 ' 



which, again, is small for M (first) and N (then) sufficiently large. This completes 
the proof. □ 



Remark 6.3. The assumption ui < has been used in (|90[) . to obtain bounds for 
Jn4>n- When the processes are stopped, as in the part of Theorem 13.21 concerning 
the case < cj < |, those estimates are essentially trivial because of the uniform 
boundedness of the stopped processes. 
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6.4 Identification of the Limiting Generator and Convergence 

In this Section we complete the proof of Theorem 13.21 The argument follows that 
of Section I5T41 so most details are omitted. 

The candidate for the limiting generator in (p2")) has been obtained in Proposition 
16.21 for the drift part, while the diffusion part comes from the term L^tfj that, by 
direct computation, is shown to be equal to 



In what follows we denote by J the generator of the diffusion process in Theorem 
13.21 The proof of convergence develops along the following steps. 
Step 1: tightness of the processes and V± . We use conditions ([50]) and 

([5Tj) . Due to (|g2")l . we are allowed to stop the processes at Tjf for some large M. 
Condition (|50l) can be obtained simultaneously for V^ 1 '^ and V^ 2 '^ by ([M]) . In 
order to establish (j5"Tj) for, e.g., V"/ 1 , we consider the function 



(i,JV) 



^(V^ N) M 2 - N) ) :=VI 

together with its perturbation 4>n as in (j5"9")) and (|FT|) . Up to o(l) terms, for stopping 
times n < T2, 

As in the proof of Proposition 16.31 we find a (possibly M-dependent) constant C 
such that the uniform bound 

i^i + "*E(^) 2 <c 



holds. This implies 



sup E 

Tl<T 2 <Tl+S 



V^ N \r 2 )-V^ N \ Tl ) 



< cs 



that, by Chebischev inequality, yields (jSTj) for Vi' N K 

Step 2: convergence to the solution of a martingale problem. Denote by [v^ 1,n \ V-[ 2 ' n ' 

a convergent subsequence of (v^ 1 '^, V± 2 ' N ^ . For a function tp : R 2 — > K of class 

C 2 and with bounded derivatives, denote by ip n its perturbation as in (|59p and (|61[) . 
Consider the martingale 

M„(*) := Vn(*) - ^n(O) - / JnMs)ds = ni J2b~ Ms)dWj(s). (93) 



It should be recalled that Vvi is a function of V^*\ V^'^, i = 1, 2, 3, 4, so when we 
write ip n (t) we mean that £ is the time at which the processes in the argument of 
tp n are evaluated. Considering that: 

• VVi — > ip as 71 ~~ ^ +°° uniformly on compact sets; 

• the processes I// 1 '™), V r 2 (l '™ ) admit a weak limit , , which is zero for 
V^' n) , i = 3,4 and V^' n \ i = 1,2,3,4, 
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it follows that the process A4 n (t) converges weakly to 



M(t) = tp (yW(t),vV\t)) - ^(^(o),^ 2 >(o)) - f Ji> (v^w.v^w 



ds. 



If we show that, for each tp with the properties specified above, M(t) is a martingale, 
then we have that the limiting processes \Vp~\t), Vi 2 \t) \ solve the martingale 
problem for J; since uniqueness holds for this martingale problem, the proof of 
Theorem I3.2I would be completed. It is therefore enough to show that M(t) is a 
martingale. Similarly to what we have done in Lemma I5.3I it suffices to show that, 
for every t > 0, 



supE 



(Mn(t)Y 



Note that 



E 



(M n (t)y 



»*£ / E 



< +00. 



d 







ds . 



Thus, it is enough to show that, for some constant C > 0, the inequality 

2" 



E 



< 



C 



n 2 



(94) 



is satisfied. 

It should be noticed that in (I92I) we gave a pointwise estimate (i.e. not in mean) 
of this sort; that, however, holds for the unperturbed function tp. In that case the 
difference between tp and its perturbation tp n was estimated by a bound of the form 
C jyP ■ But now we are not stopping the process anymore, so a little more care is 
needed. We recall that 

tp n = tp + n~itpx. 

Given the bound in in order to obtain (|M|) it is enough to show that 

2" 



E 



C 
< — . 

n 



(95) 



As seen in Section 16.21 tpi is a linear combination of terms of the form 

F := ditp (Jv^dpnijvf^dpn^j J v[ l) dp n J v { 2 h) dp n , 
i = 1, 2, I, h = 1, 2, 3, 4. So it is enough to consider one of such terms. We have 



dF 

r: :i 



dx i n 4 



( / v^dpn, I vf'dp n ) jj^v\ lJ (xj,Tij) / v\ J ' dp n j v^'dp n 



,(2) 



9 ..(1)/ 



..CO. 



,00, 



+< 9 2, l V> ^ v[ 1} dp n , J v[ 2) dp n ^j ^[^{x^rij) J v[ 

+ dttp (^J v^dpn, J v^dpr^j -^v[ l \xj,rij) J vi, h) dp. 



+ditp \J v^dpn, J vf'dp n J ■^v^'{xj,rij) J v[ l 'dp„ 

Consider the first of the summands above, the others can be dealt with similarly. 
The factor 



,( 2 ) 



,00, 



d 'ii^ [ I v[ 1] dp n , j V^'dpn ) j^vY'ixj,^) 
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is uniformly bounded. Also the term 



4 / V^dP; 
714 J 



•n 



is uniformly bounded. The last factor, J dp n , is clearly bounded in absolute 
value by ||/5 n || r) defined in (|66[) . Estimating similarly all terms, one sees that 



for some constants C, C", where we have used (|5Tj) . This establishes (|9"5"j) , and thus 
completes the proof of Theorem 13.21 
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